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Abstract 

We report measurements of radiative B decays with if 777 final states, using a data 
sample of 253 fb -1 recorded at the T(4S) resonance with the Belle detector at 
the KEKB e + e~ storage ring. We observe B + — > K + r]^ for the first time with a 
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branching fraction of (8.4 ± 1.5(stat) ^(syst)) x 10 ~ 6 for M K V < 2.4 GeV/c 2 , and 
find evidence of B° — > K i]^. We also search for i? — > 1£|(1780) 7 . 

Key words: radiative B decay 
PACS: 13.20.He, 14.40.Nd 



Radiative B decays, which proceed mainly through the b — > 57 process 2 , 
have played an important role in a search for physics beyond the Standard 
Model (SM). Although the inclusive branching fraction has been measured to 
be (3.3 ±0.4) x 10 -4 [1], we know little about its constituents. So far, measured 
exclusive final states such as If* (892)7 [2,3], K* (1430)7 [2,4], Kirir-f [4] and 
K<fyy [5] only explain one third of the inclusive rate. Detailed knowledge of 
exclusive final states reduces the theoretical uncertainty in the measurement 
of the inclusive B — > X s 7 branching fraction using the pseudo-reconstruction 
technique, as well as in the measurement of B — > X s l + l~ [6]. In this anal- 
ysis, the decay mode B — > Kr/j is studied for the first time. In addition to 
improving the understanding of b — > 57 final states, B° — > Kgr/j can be used 
to study time-dependent CP asymmetry [7], which is sensitive to physics be- 
yond the SM. The mode B — > Krj'-f can also be used to search for radiative B 
decays through possible Krj resonances, e.g., i^g(1780) observed by the LASS 
experiment [8]. 

The analysis is based on 253 fb" 1 of data taken at the T(AS) resonance 
(on-resonance) and 28 fb _1 at an energy 60 MeV below the resonance (off- 
resonance), which were recorded by the Belle detector [9] at the KEKB asym- 
metric e + e~ collider (3.5 GeV on 8 GeV) [10]. The on-resonance data corre- 
sponds to 275 million BB events. The Belle detector is comprised of a silicon 
vertex detector, a 50-layer central drift chamber (CDC), an array of aerogel 
Cherenkov counters (ACC), time-of-flight scintillation counters (TOF) and an 
electromagnetic calorimeter of CsI(Tl) crystals (ECL) located inside a super- 
conducting solenoid coil that provides a 1.5 T magnetic field. An instrumented 
iron flux-return for K^f/j, detection is located outside the coil. Two different 
inner detector configurations were used. For the first sample of 152 million 
BB pairs, a 2.0 cm radius beampipe and a 3-layer silicon vertex detector were 
used; for the latter 123 million BB pairs, a 1.5 cm radius beampipe, a 4-layer 
silicon detector and a small-cell inner drift chamber were used [11]. 

We reconstruct B + — > K + r] r ) and B° — > Kgi]'-/ via i] — > 77 and 77 — > tt + 7t~7t°. 
All charged tracks used in the reconstruction (except charged pions from 
Kg) are required to have a center-of-mass (CM) momentum greater than 

1 on leave from Nova Gorica Polytechnic, Nova Gorica, Slovenia 

2 Throughout this paper, the inclusion of the charge conjugate mode is implied 
unless otherwise stated. 
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100 MeV/c and to have an impact parameter within ±5 cm of the interaction 
point along the positron beam axis and within 0.5 cm in the transverse plane. 
In order to identify kaon and pion candidates, we use a likelihood ratio based 
on the light yield in the ACC, TOF information and specific ionization mea- 
surements in the CDC. For the requirement applied on the likelihood ratio, 
we obtain an efficiency (pion misidentification probability) of 90% (10%) for 
charged kaon candidates, and an efficiency (kaon misidentification probability) 
of 98% (10%) for charged pion candidates. Tracks identified as electrons are 
excluded. 

K s candidates are formed from it + it~ combinations with invariant mass within 
8 MeV/c 2 (~ 2a) of the nominal K° s mass. The two pions are required to have 
a common vertex displaced from the interaction point. The Kg momentum 
direction is required to be consistent with the K s flight direction. Neutral 
pion candidates are formed from pairs of photons that have an invariant mass 
within 16 MeV/c 2 (~ 3a) of the nominal n° mass and a momentum greater 
than 100 MeV/c in the CM frame. Each photon is required to have an energy 
greater than 50 MeV in the laboratory frame. A mass-constrained fit is then 
performed to obtain the 7T° momentum. 

For i] — > 77 reconstruction, we require that the invariant mass of the two 
photons satisfy 0.515 GeV/c 2 < M 77 < 0.570 GeV/c 2 and that each photon 
have an energy greater than 50 MeV in the laboratory frame. We also require 
| cos#{J el | < 0.9, where #{J el is the angle between the photon momentum and 
i] boost direction from the laboratory frame in the rj rest frame. A mass- 
constrained fit is then performed to obtain the 77 momentum. For 77 — > 7r + 7r~7r°, 
we apply a requirement on the three-pion invariant mass, 0.532 GeV/c 2 < 
M n+n - n o < 0.562 GeV/c 2 . 

We reconstruct B meson candidates from an 77, a charged or neutral kaon and 
the highest energy photon within the acceptance of the barrel ECL (33° < 
# 7 < 128°, where # 7 is the polar angle of the photon with respect to the 
electron beam in the laboratory frame). Here, the invariant mass of the Kr] 
system is required to be less than 2.4 GeV/c 2 . This selection corresponds to 
E^ > 2.1 GeV, where is the photon energy in the B rest frame, and 

includes 84% of events from the b — > S7 process. The highest energy photon 
candidate is required to be consistent with an isolated electromagnetic shower, 
i.e., 95% of the energy in an array of 5 x 5 crystals should be concentrated 
in an array of 3 x 3 crystals and no charged tracks should be associated with 
it. In order to reduce the background from decays of 7T° and 77 mesons, we 
combine the photon candidate with each of the other photons that have CM 
energy greater than 30 MeV (200 MeV) in the event and reject the event if the 
invariant mass of any pair is within 18 MeV/c 2 (32 MeV/c 2 ) of the nominal 
7T° (77) mass. This condition is referred to as the TT°/r] veto. 
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We use two independent kinematic variables for the B reconstruction: the 
beam-energy constrained mass M bc = \J (E^ eam / c 2 ) 2 — (\pk v + P^\/ c ) 2 an d 
AE = E* Kr) + E* - £ beam , where E£ cam is the beam energy, and p*, E*, 
Px v , E* Krj are the momenta and energies of the photon and the Krj system, 
respectively, calculated in the CM frame. In the M bc calculation, the photon 
momentum is rescaled so that \p*\ = (-E bcam — E* Kr] )/c is satisfied. We require 
M bc > 5.2 GeV/c 2 and -150 MeV < AE < 80 MeV. We define the B signal 
region to be M bc > 5.27 GeV/c 2 . In the case that multiple candidates are 
found in the same event, we take the candidate that has the 77 mass closest 
to the nominal mass 3 after applying the background suppression described 
later. 

The largest source of background originates from continuum e + e" — > qq 
(q = u, d, s, c) production including contributions from initial state radiation 
(e + e~ — > qqj). In order to suppress this background, we use the likelihood 
ratio (LR) described in Ref. [4], which utilizes the information from a Fisher 
discriminant [12] formed from six modified Fox- Wolfram moments [13] and 
the cosine of the angle between the B meson flight direction and the beam 
axis. The LR requirement retains 44% of the signal, while rejecting 98% of 
the continuum background. 

In order to extract the signal yield, we perform a binned likelihood fit to the 
Mb c distribution. The M bc distribution of the signal component is modeled 
by a Crystal Ball line shape [14], with the parameters determined from the 
signal Monte Carlo (MC) and calibrated using control samples of B + — > D°(^ 
K + n-n°)n+ and B° -> K%n-n°)n + decays. The M bc distribution of 

the continuum background is modeled by an ARGUS function [15] whose 
shape is determined from the off-resonance data. Here, the LR requirement is 
not applied to the off-resonance data in order to compensate for the limited 
amount of data in that sample. The possible bias due to this is taken as 
systematic error on the fitted yield. Background from hadronic B decays is 
divided into two components, which we refer to as generic BB background and 
rare B background in this paper. The former comprises B decays through b — > 
c transitions including color-suppressed B decays such as B° — > D tc°, and the 
latter covers charmless B decays. Each of them is modeled by another ARGUS 
function. The shape of these distributions is determined using corresponding 
MC samples. In order to study the contamination from other b — > 57 decays, 
we examine a B — > K* (892)7 MC sample and an inclusive b — > 57 MC sample 
that is modeled as an equal mixture of sd and su quark pairs and is hadronized 
using JETSET [16], where the X s mass spectrum is fitted to the model of 
Kagan and Neubert [17]. We find that feed-down from other b — > 57 decays 
is small, but not negligible, and model its M bc distribution with an ARGUS 



In case multiple candidates share such an r] candidate, the candidate with the 
smallest |A£1 is chosen. 
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Fig. 1. M hc distributions for (a) B + — > K + r/j, (b) 5° — > KgTjj. Fit results are 
overlaid. 



function. 



Figure 1 shows the M^ c distributions for B + — > K + r] r y and 5° — > Kgrj^y, re- 
spectively. These distributions, as well as the distribution for the combined 
mode, are fitted to the sum of signal, continuum, generic BB, rare £? back- 
ground and b — > S7 feed-down components. In the fit, the normalization of 
generic 55, rare 5 and 6 — > 57 are fixed according to the luminosity and 
b — > 57 branching fraction, while the normalization of the continuum compo- 
nent is allowed to float. We find signal yields of 81 ± 14, 20.9 tH and 102 ± 16 
events with statistical significances of 7.1a, 3.7a and 8.1a", for the charged, 
neutral and combined modes, respectively. Here, the significance is defined 
as ^—2 ln(£(0)/£ max ), where £ max and £(0) are the maximum values of the 
likelihood when the signal yield is left free or fixed to zero, respectively. 

Figure 2 shows the 77 and 7r + 7r~7r° invariant mass distributions for events 
inside the B signal region. Here, we do not apply the best candidate selec- 
tion. We observe clear peaks at the nominal r\ mass. The Kr\ invariant mass 
distribution for events inside the B signal region is shown in Fig. 3. Here, the 
background distributions are obtained from the off-resonance data without the 
LR requirement or from the corresponding MC samples, and are normalized 
using the fit result. We find that the signal is concentrated between 1.3 GeV/c 2 
and 1.9 GeV/c 2 and is falling above 1.9 GeV/c 2 . Therefore, our requirement 
Mkt] < 2.4 GeV/c 2 is expected to include most of the B — ► Kr]^ signal. We 
do not see any clear resonant structure in the M Krj distribution. 

The systematic error on the signal yield due to the fitting procedure is esti- 
mated by varying the value of each fixed parameter by ±la and extracting the 
new signal yield for each case. The difference between the background shape 
for the continuum MC with and without the LR requirement is taken as an 
additional error to the continuum background shape. We set the normaliza- 
tion of either the generic BB or rare B backgrounds to zero and to twice its 
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Fig. 2. 77 invariant mass distributions for (a) 77 — > 77 and (b) 77 — > 7r + 7r~7r° inside the 
signal region for combined B — > If 777. Dashed lines show the selection applied in 
the analysis. 
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Fig. 3. X77 invariant mass distribution for events inside the B signal region for 
combined B — > if 777. 

nominal value to account for its uncertainty. The changes of the yields for 
each procedure are added in quadrature, and are regarded as the systematic 
error on the signal yield. We also calculate a statistical significance for each 
case, and regard the smallest value as the significance including the systematic 
error. The result is listed in Table I. 




The signal reconstruction efficiency is estimated using the MC simulation and 
is corrected for discrepancies between data and MC using control samples. The 
signal MC has uniform Kr\ invariant mass and cos #hei distributions, where #h e i 
is the decay helicity angle between the kaon momentum and opposite to B 
momentum in the Kt] rest frame. We find that the efficiency is almost indepen- 
dent of the Kr] invariant mass and cos#hei- Table 1 shows the signal efficiencies 
and the branching fractions for each B — > ^777 mode. Here, we assume an 
equal production rate for B°B and B + B~ . The error on the branching frac- 
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Table 1 

Measured signal yields, efficiencies, branching fractions (B) and significances includ- 
ing systematic error (<S) for B — ► K-qj. The first and second errors are statistical 
and systematic, respectively. Efficiencies include the sub-decay branching fractions. 





Mode 


Yield 


Efficiency (%) 


£ (xl(T 6 ) 


S 


B+ 
B° 
B - 


— > K + r]^ 
-» K°r?7 
-> K7/7 


81 ±14^° 

on q +7.3 +4.2 
^ U - y -6.5 -3.2 

102 ± 16 tf 


3.50 ±0.27 
0.87 ±0.08 
4.37 ±0.31 


8.4±1.5±S;1 

o 7 +3.1 +1.9 
°- ' -2.7 -1.6 

8.5 ±1.3^| 


6.8 
3.4 
7.7 



tion includes the following systematic uncertainties: photon detection (2.8%), 
tracking (1.0% to 1.2% per track), kaon identification (0.8%), pion identifica- 
tion (0.5% per pion), Kg detection (4.5%), 7r° detection (1.5%), i] detection 
in rj — > 77 mode (2.0%), 7i /r] veto and LR (5.9% and 4.4% for charged and 
neutral modes, respectively), possible Krj mass dependence of the efficiency 
(2.1% and 4.4% for charged and neutral modes, respectively), possible cos^hei 
dependence of the efficiency (2.5% and 3.4% for charged and neutral modes, re- 
spectively), uncertainty in the f] branching fraction (0.7% for r] — > 77 and 1.8% 
for rj — > 7T + 7r~7r ), and uncertainty in the number of BB events (1.1%). The 
systematic errors from the 7i°/r] veto and LR requirement are estimated using 
control samples of B + -> D°(^ K+n-n°)n + and B° -> £>-(-> ^°7r-7r°)7r + 
decays, treating the primary pion as a high energy photon. 

We search for the decay B — > (1780)7 by applying the additional re- 
quirements 1.60 GeV/c 2 < M Kv < 1.95 GeV/c 2 and | cos 6" he i I < 0.2 or 
I cos#hei| > 0.7. The expected cos#hci distribution for B — > ^(1780)7 is pro- 
portional to 1 — 11 cos 2 #hei + 35 cos 4 #hei — 25 cos 6 #hci- The fits to the Mb c dis- 
tributions yield 4.4 t 2 2 % 0.2 t\l and 5.2 ±f;| events for the charged, 
neutral and combined modes, respectively. Here and in the following, we quote 
statistical and systematic errors in the first and second position. The M bc dis- 
tribution and fit result for the combined mode is shown in Fig. 4. We provide 
only upper limits due to our inability to distinguish B — > (1780)7 from non- 
resonant decays. The 90% confidence level upper limit N is calculated from the 
relation C(n)dn = 0.9 J °° C(n)dn, where C(n) is the maximum likelihood 
in the M bc fit with the signal yield fixed at n. In order to include the systematic 
errors from the fitting procedure in the upper limit for the yield, the positive 
systematic error is added to N. The obtained yield upper limits, efficiencies 
and products of branching fractions B(B — > (1780)7) x B(K^(1780) — > Kr\) 
are listed in Table 2. Here, the number of BB events and the reconstruction 
efficiency are lowered by la when we calculate the upper limit for the branch- 
ing fractions. If we assume B(K*(1780) -> Kr}) = (11 ±I)% [18], the 90% 
confidence level limits correspond to B — > ii^ (1780)7 branching fractions of 
3.9 x 10~ 5 , 8.3 x 10~ 5 and 3.7 x 10~ 5 , respectively for charged, neutral and 
combined modes, which substantially improve the limits set by the ARGUS 
collaboration [19]. 



9 



"si 20- _r~L data 

| total 

w — continuum background 

|T - background from B 

§j I b->sy feed-down 

13 " 

10 - n 




M bc (GeV/c 2 ) 

Fig. 4. M bc distribution for combined B — > ^777 with the B — > Kg (1780)7 selection. 
Fit results are overlaid. 

Table 2 

Measured signal yields, efficiencies and products of branching fractions of B — > 
#1(1780)7 and X|(1780) Kr] (B x B{K% -» K77)). Efficiencies include the sub- 
decay branching fractions of 77 and if , but not of if|(1780). Upper limits are 
calculated at the 90% confidence level and include systematics. 





Mode 


Yield 


Efficiency (%) 


B x B{Kl^Kr]) (xlO" 6 ) 


B+ 


^if 3 *(1780)+ 7 


< 15.0 


2.03 ±0.16 


< 2.9 


B°- 


-> K|(1780)°7 


< 7.5 


0.48 ± 0.05 


< 6.4 


B - 


^1(1780)7 


< 17.7 


2.51 ±0.18 


< 2.8 



Some extensions of the SM predict a large CP asymmetry in the b — > S7 pro- 
cess [20]. We measure the partial rate asymmetry = (1/(1 — 2u>))(iV_ — 
iV + ) / (AL±iV + ) for B + — > K + i]^f, where is the signal yield for B T — > K T r]^f 
and u> is the probability that a signal event is reconstructed with the wrong 
kaon (and hence B) charge. This probability is found to be less than 1% 
in our signal MC sample, and hence we ignore its negligible effect on A C p. 
N T is obtained by fitting separately the M^ c distributions for the negatively 
and positively charged modes shown in Fig. 5. We find iV_ = 34.0 l;^ and 
N + = 46.7 tg° 8 5 . The systematic error on Acp consists of the following con- 
tributions. The error from the fitting procedure is estimated to be 0.045 by 
varying each fixed parameter one by one, and extracting Acp for each pro- 
cedure, in the same way as before. Here, we assume no asymmetry for the 
generic BB background, but allow 100% asymmetry for the rare B and 6% 
asymmetry for b — > 37 [21]. The error from the overall detector bias is studied 
with the B° — > D~(K~it + -k )-k + control sample and is found to be 0.035. By 
adding these errors and the possible asymmetry in kaon identification (0.014) 
in quadrature, we obtain A C p = —0.16 ± 0.09 ± 0.06. 

In conclusion, we observe the decay mode B + K + i] , ~f and find the first 
evidence of B° K°rj^. The branching fraction and partial rate asymmetry 
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Fig. 5. Mbc distributions for (a) negative charged B 
B + — > K + 7]^j. Fit results are overlaid. 

of B + -> i^+r?^ are measured to be (8.4±1.5 ^j) xlO" 6 and -0.16±0.09±0.06 
for Mkti < 2.4 GeV/c 2 . The branching fraction of B° — > K°i] / -f is measured 



to be (8.7 



+3.1 +1.9\ 
-2.7 -1.6J 



x 10 . We also search for B — > (1780)7, but find no 
evidence. Although the signal yield for B° — ■> Kgi]'-/ is small, this mode can be 
used in the near future to study time-dependent CP asymmetries in radiative 
B decays and to search for new physics. 
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